JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOL10.1029/, 



Kinetic-scale magnetic turbulence and finite Larmor 
radius effects at Mercury 

V. M. Uritsky\ J. A. Slavin^ G. V. Khazanov^ E. F. Donovan\ S. A. 

2 3 3 

Boardsen , B. J. Anderson , and H. Korth 



V. M. Uritsky, Physics and Astronomy Department, University of Calgary, SB605, 2500 Uni- 
versity Drive NW, Calgary, AB T3A0P4, Canada. (uritsky@ucalgary.ca) 

^University of Calgary, Calgary, AB, 
Canada 

^NASA Goddard Space Flight Center, 
Greenbelt, MD, USA 

^John Hopkins University Applied 
Physics Laboratory, Laurel, MD, USA 



DRAFT 



January 19, 2013, 5:52pm 



DRAFT 



X - 2 URITSKY ET AL.: KINETIC-SCALE TURBULENCE AT MERCURY 

3 Abstract. We use a nonstationary generalization of the higher-order struc- 

4 ture function technique to investigate statistical properties of the magnetic 

5 field fluctuations recorded by MESSENGER spacecraft during its first flyby 

6 (01/14/2008) through the near Mercury's space environment, with the em- 

7 phasis on key boundary regions participating in the solar wind - magneto- 

8 sphere interaction. Our analysis shows, for the first time, that kinetic-scale 

9 fluctuations play a significant role in the Mercury's magnetosphere up to the 

10 largest resolvable time scale (~ 20 s) imposed by the signal nonstationar- 

11 ity, suggesting that turbulence at this planet is largely controlled by finite 

12 Larmor radius effects. In particular, we report the presence of a highly tur- 

13 bulent and extended foreshock system filled with packets of ULF oscillations, 

14 broad-band intermittent fluctuations in the magnetosheath, ion-kinetic tur- 

15 bulence in the central plasma sheet of Mercury's magnetotail, and kinetic- 

16 scale fluctuations in the inner current sheet encountered at the outbound (dawn- 

17 side) magnetopause. Overall, our measurements indicate that the Hermean 

18 magnetosphere, as well as the surrounding region, are strongly affected by 

19 non-MHD effects introduced by finite sizes of cyclotron orbits of the consti- 

20 tuting ion species. Physical mechanisms of these effects and their potentially 

21 critical impact on the structure and dynamics of Mercury's magnetic field 

22 remain to be understood. 
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1. Introduction 

23 Dynamic variability of Mercury's magnetosphere has been intensively studied in the con- 

24 text of tail and magnetopause reconnection, magnetic flux transport, ULF waves and oscil- 

25 lations, and other phenomena (see e.g. Anderson et al. [2008]; Slavin et al. [2008, 2009a]; 

26 Boardsen et al. [2009a]; Sundherg et al. [2010]). However, little is known about mag- 

27 netic turbulence in the Hermean plasma environment. Korth et al. [2010] investigated 

28 turbulence in the unperturbed solar wind observed by MESSENGER at the hehocentric 

29 distances of Mercury's orbit, but they did not address magnetic fluctuations formed in 

30 the vicinity of the planet. Considering the predicted significance of finite Larmor ra- 

31 dius (FLR) effects and the abundance of plasma instabilities in Mercury's magnetosphere 

32 [Glassmeier and Espley, 2006; Blomberg et al, 2007; Travnicek et al., 2009], one can ex- 

33 pect the magnetic fluctuations at this planet to be heavily affected by plasma kinetics. 

34 This prediction, however, has not been tested on empirical data until now. 

35 Turbulence in plasma is a fundamental physical phenomenon in its own right (see, for 

36 instance, Pouquet [1978]; Politano and Pouquet [1995]; Khazanov et al. [1996]; Robinson 

37 [1997]; Biskamp [2003]; Mininni and Pouquet [2007]; Singh et al. [2007]; Schekochihin 

38 et al. [2009]). Large-scale stochastic plasma motions are controlled by the magnetohydro- 

39 dynamic (MHD) energy cascade involving Alfvenic wave packets of various sizes which 

40 tend to violate statistical laws derived for nonmagnetized fluids by exhibiting strong spa- 

41 tial anisotropy [Schekochihin et al., 2007]. In the resistive MHD approximation, the width 

42 of the (inertial) range of scales defined by this regime is reflected by the magnetic Reynolds 

43 number, with the main dissipation taking place at the distances shorter than the Taylor 
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44 microscale [Uritsky et ai, 2010a]. The situation is substantially more complicated in col- 

45 lisionless plasmas with vanishing resistivity where the upper cutoff of the inertial range 

46 in the wave-number space is created by ion kinetics. The latter usually generates a new 

47 cascade with distinct physical and statistical properties, transferring the energy to even 

48 smaller scales. Yet another type of turbulence is found at the electron scales as demon- 

49 strated recently for the solar wind [Sahraoui et al, 2009]. Understanding these effects is 

50 an important problem with significant theoretical implications, including the fundamental 

51 mechanisms of magnetic reconnection, particle acceleration and transport in space plas- 

52 mas [Chang, 1999; Lazarian and Vishniac, 1999; Antonova, 2002; Borovsky and Funsten, 

53 2003; Untsky et a/., 2001; Pulkkinen et a/., 2006; Uritsky et a/., 2008; Servidio et a/., 2009; 

54 Eastwood et al, 2009; Klimas et al, 2010]. 

55 In this study, we pursue a more practical goal by applying turbulent analysis tools as 

56 a means of characterizing ion kinetic scales in different plasma structures surrounding 

57 Mercury. Our methodology is based on the existence of scaling crossover separating 

58 MHD and ion kinetic regimes of magnetic fluctuations. By identifying this crossover in 

59 the temporal domain and mapping the results to the wave-number space using predicted 

60 values of flow velocity, we evaluate the ion gyro radius and temperature in several locations 

61 of the Hermean magnetosphere, and compare these measurements with earlier theoretical 

62 estimates. We also show that scaling regimes of magnetic fluctuations vary greatly in 

63 the Mercury's foreshock, magnetosheath, and the magnetosphere, and that they involve 

64 contributions from a variety of non-random processes and structures, including boundary 

65 layers, rotating flows, and transient ULF activity. Overall, our analysis suggests that ion 

66 kinetic turbulence is present in all Hermean plasma structures, and is the leading source of 
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67 stochastic variability inside the magnetopause up to the largest resolvable scales imposed 

68 by a data nonstationarity. 

69 The paper has the following structure. The next section briefly summarizes properties of 

70 magnetic fluctuations essential for this research, explains the link between the Fourier and 

71 structure function analyses, and discusses several classes of turbulent cascades. Section 3 

72 describes the methods and the data used in this study. One of the mathematical tools, 

73 the continuous structure function scalogram, is introduced for the first time and to our 

74 knowledge has not been used in space or turbulence studies before. In section 4, we present 

75 the results of the analysis of magnetic fluctuations recorded during MESSENGER'S flrst 

76 flyby. We investigate nonstationary scaling structure of the these fluctuations and provide 

77 comparative description of turbulent regimes in several key plasma regions visited by the 

78 spacecraft. Finally, a summary of plasma parameters assessed using the measured ion 

79 crossover scales is presented and discussed in the context of previous investigations. 

2. Brief theoretical background 

80 The autocorrelation properties of turbulent fluids are commonly described in frames of 

81 two complementary statistical formalisms, the Fourier analysis and the structure function 

82 approach [Politano and Pouquet [1995]). 

The time-domain higher-order structure function (SF) is deflned as 

s,{t) = , (1) 

83 in which 5Bt are the increments of the studied turbulent fleld B measured at time lag r, 

84 (•) denotes averaging over all pairs of points separated by this lag, and q is the order. The 

85 SF exponents Q estimated from the scaling ansatz Sq{T) oc r^«, along with the spectral 
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86 exponent (3 describing the power-law decay of the wave-number Fourier power spectrum 

87 P{k) oc k"^, provide a detailed description of the turbulent regime under study. The 

88 second-order SF 5*2 (r) plays a special role in statistical mechanics of turbulent media as 

89 a proxi to the band- integrated spectrum [Biskamp, 2003], yielding (2 — P — ^ under the 

90 assumption of linear space-time coupling as will be discussed later. 

91 Fig. 1 illustrates the overall shape and mutual relationship between the wave-number 

92 Fourier power spectrum and the temporal SF for a typical turbulent environment observed 

93 by a spacecraft. Both statistical descriptions reveal three fundamental scaling regions 

94 labeled I through III in the Figure. Each region is characterized by its own set of power 

95 laws with distinct values of spectral and SF exponents. 

96 The large-scale region I is dominated by MHD and hydrodynamic energy cascades 

97 represented by "fluid" exponent values /3 ~ 1.5 — 2.0 and ^2 ~ 0.5 — 1.0. The frequently 

98 observed ^ = 5/3 (C2 = 2/3) corresponds to the Kolmogorov scahng of the Alfven-wave 

99 energy spectrum of perpendicular wave modes indicative of a fully developed turbulent 

100 state [Biskamp, 2003]. The same 5/3 law describes the hierarchy of isotropic turbulent 

101 eddies in non-magnetic fluids [Kolmogorov, 1941]. As the wave number grows tending 

102 to the inverse ion gyro radius p^^, plasma kinetics becomes increasingly important. For 

103 kpi < 1 (or equivalently for r > r^, where Tj is the position of the ion crossover in 

104 the time-lag space), the kinetic effects can be treated by an extended MHD approach 

105 with kinetically-calculated anisotropic pressure tensor [Schekochihin et al, 2007]. In this 

106 sub-kinetic regime, the energy cascade continues to be supported by counter-propagating 

107 Alfven wave packets, and the scaling exponents are approximately the same as in the 

108 usual MHD regime. 



DRAFT 



January 19, 2013, 5:52pm 



DRAFT 



URITSKY ET AL.: KINETIC-SCALE TURBULENCE AT MERCURY X - 7 

109 At the ion crossover scales kpi ~ 1 and t Ti, the kinetic MHD approximation breaks 

no down since the Alfvenic fluctuations are no longer decoupled from the kinetic compo- 

111 nent of the turbulence represented by density and magnetic-field strength fluctuations 

112 [Schekochihin et al, 2007]. The resulting scaling regime II, which will be referred to 

113 as the ion-kinetic regime throughout this paper, is characterized by the "ion" values of 

114 spectral and SF exponents which are larger than the their fluid counterparts, leading to 

115 steeper log-log slopes of P{k) and Sg{T). The micro turbulence theories developed for this 

116 range of scales predict /3 !v 2.3 — 2.5 (or ^2 ~ 1-3 — 1.5), depending on the underlying 

117 dispersive wave mode (usually kinetic Alfven waves (KAW) or whistler branches with 
lis secondary lower hybrid activity), and the turbulence type (i.e. a weak or strong), see 

119 e.g. Yordanova et al. [2008]; Eastwood et al. [2009]; Sahraoui et al. [2009]. Compressional 

120 corrections tend to increase the exponents [Alexandrova et al., 2008] making them deviate 

121 further from fluid values. 

122 As k reaches the inverse electron gyro radius (mapped to the time scale Te), both 

123 electrons and ions become demagnetized leading to a new scaling regime III with even 

124 higher /3 and ( values. Physical mechanism of the electron cascade is currently not well 

125 understood, with obhque KAW modes being a candidate explanation for the cross-scale 

126 coupling in this regime [Sahraoui et al., 2010]. 

127 Due to the coUisionless dissipation at ion and electron gyroscales, only a certain fraction 

128 of the turbulent power in regions II and III arriving there from fluid scales is converted 

129 into turbulent cascades, while the rest is subject to Landau damping and other types 

130 of wave-particle interaction [Khazanov, 2010], steepening the apparent spectral and SF 
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131 slopes and making the kinetic-scale turbulence even more distinguishable from the fluid 

132 regime. 

3. Data and Methods 

133 Keeping this theoretical framework in mind, we investigated time series of magnetic 

134 field variations recorded by the MAG instrument [Anderson et ai, 2007] onboard the 

135 MESSENGER spacecraft during its first flyby near Mercury (closest approach at 19:04:39 

136 01/14/2008, sampling frequency 20 Hz). The northward direction of the interplanetary 

137 magnetic field during this flyby provided ideal conditions for studying intrinsic properties 

138 of Mercury's magnetic fluctuations not distorted by a pronounced susbtorm activity. As 

139 we show below, MESSENGER MAG is able to resolve the ion kinetic scales in all of the 

140 Hermean plasma structures visited during the first flyby. The sampling frequency of MAG 

141 is also sufficient for identifying these scales in the surrounding solar wind as demonstrated 

142 by KoHh et al. [2010]. 

143 Fig. 2 shows the trajectory of the spacecraft relative to the average positions of the 

144 Mercury's bow shock and magnetopause [Slavin et al., 2009a]. Using the nonstationary 

145 data analysis tools described below, we studied the evolution of magnetic turbulence 

146 during the entire flyby, and performed a more focused investigation of selected plasma 

147 regions discussed in previous publications [Slavin et al., 2008, 2009a; Boardsen et al., 

148 2009a; Sundherg et al., 2010]. These regions represent important boundary layers and 

149 processes forming the response of Mercury's magnetosphere to the solar wind driver. 

150 The following regional identifiers are used throughout the paper: SWl - unperturbed 

151 solar wind at the dusk side; FSl - outermost dusk-side foreshock region; FS2 - innermost 

152 foreshock region near the dusk bow shock; MSI - outermost magnetosheath at the dusk 
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153 flank; FTE - one-minute interval involving a flux transfer event [Slavin et ai, 2008]; MS2 

154 - innermost magnetosheath contacting the dusk magnetopause; KH - Kelvin - Helmholtz 

155 activity inside the dusk magnetopause [Slavin et ai, 2008; Sundberg et ai, 2010]; CCS - 

156 cross-tail current sheet; DD - first diamagnetic decrease encountered in the inner mag- 

157 netosphere; IBL - ion boundary layer adjacent to the dawn magnetopause [Slavin et al., 

158 2008; Anderson et al, 2011]; MS3 - innermost magnetosheath observed immediately after 

159 exiting dawn magnetopause; MS4 - outermost magnetosheath before crossing the dawn- 

160 side bow shock; FS3 - innermost foreshock region on the dawn side ; FS4 - outbound 

161 foreshock adjacent to the unperturbed solar wind; SW2 - solar wind observed at the end 

162 of the flyby. Timing information for each region is provided in Table 1. For the reader's 

163 convenience, the regions are also marked by color-coded bars in Figs. 2 and 3. 

164 To identify the ion crossover scales and other characteristic features of Mercury's mag- 

165 netic turbulence, we used the method of higher order structure functions generalized for 

166 the case of strongly nonstationary signals. 

167 Our choice of the SF-based approach to the magnetic turbulence at Mercury, as opposed 

168 to a somewhat more popular Fourier analysis, is motivated by two factors. First, SF 

169 analysis is a powerful tool of differentiating between random and deterministic components 

170 of multiscale variability. While the second-order SF contains essentially the same scaling 

171 information as the power spectrum (see Section 2), the other SF orders {q ^ 2) provide 

172 additional important clues on the structure of the studied signal; in particular, they enable 

173 identiflcation of non-random transient disturbances mixed with stochastic noise. A purely 

174 random signal is described by the condition dC,q/dq > 0, with the functional dependence 

175 of ( on the order q quantifying the stochastic intermittency ( "spikiness" ) of the data 
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176 [Politano and Pouquet, 1995]. A set of "flat" SFs with Q ^ Wq reveals the presence 

177 of singular features such as discontinuities and shocks. The inverted hierarchy of SF 

178 exponents {d(g/dq < 0) is indicative of a (quasi-) periodic wave oscillation embedded in 

179 a stochastic background, with the largest scale satisfying this criterion giving the period 

180 of the oscillation. 

181 Secondly, the SF analysis is generally more robust when applied to nonstationary signals, 

182 as well as when the data amount is scarce. This advantage is crucial because MESSEN- 

183 GER's magnetic measurements contain strong trends reflecting spatial inhomogeneity of 

184 the traversed plasma structures. Because of these trends, direct calculation of spectral 

185 power from a Fourier transform can be quite inaccurate, especially at the low frequencies 

186 comparable with the inverse time scale of the trends. Since the trends are nonlinear, de- 

187 trending the data introduces uncontrolled spectral errors and does not resolve the problem. 

188 SF analysis is much less sensitive to such effects and is statistically more stable when ap- 

189 plied to short data sets, the properties that are particularly useful for a windowed analysis 

190 of fly by time series. 

191 The presence of scaling crossovers is usually evident in both the wave-number and the 

192 time-lag representations as illustrated by Fig. 1. However, the mapping between the 

193 crossover scales as seen in the k and r domains depends on the state of the plasma and is 

194 not always straightforward. In the simplest case, when the bulk flow velocity vq is much 

195 higher than the characteristic propagation speed of the wave modes underlying turbulent 

196 motion, the Taylor "frozen-in flow" approximation cu = kvo can be applied, which yields 

197 k — 2t:/vqt. By applying the ion-crossover condition kpi ~ 1, we can therefore evaluate 
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198 



the ion gyro radius and the ion temperature (in energy units): 



Pi ~ VoTi/2TT, 



(2) 



Ti fa mi{voTi/Tci)'^, 



(3) 



199 in which Tj is the ion crossover scale obtained from the temporal SF analysis, Cj and 

200 nii are respectively the charge and the mass of the ions, Td — 2T:mi/eiB is the local 

201 gyroperiod, and B is the local magnetic field. The Taylor assumption used in these 

202 relations is approximately valid for the solar wind, the magnetosheath, the magnetopause 

203 boundary, and for the magnetotail plasma sheet [Matthaeus et al., 2005; Alexandrova 

204 et al, 2008; Yordanova et al, 2008; Vows et al, 2006]. For other magnetospheric regions, 

205 this assumption can be inapplicable and the space-time coupling far from trivial. 

206 To deal with signal nonstationarity, we used the sliding window technique. The time 

207 series under investigation was segmented into a sequence of overlapping intervals of a 

208 fixed width A representing an empirical compromise between the nonstationarity and the 

209 intrinsic variability of the data, shifted by a constant shift A/2. For each window position, 

210 we computed a set of SFs according to eq. 1 given in Section 2, with q e {1, 2, 3, 4} and 



The time-dependent shape of the resulting two-dimensional windowed structure function 
Sg{T, t), with t being the running time variable given by the central position of the sliding 
window, was represented in two different formats: (1) by the time series (q{t) of scaling 
exponents estimated over several selected r ranges, and (2) as the continuous time - period 
scalogram Cqi'^^'t) enabling classification of turbulent regimes across the entire range of 



211 



T < A/2. 
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available time scales: 



d log 




91 



d log T 



(4) 



212 Here, B{t') = B{t') — 4>{t, t', A) is the locally detrended magnetic signal, is the quadratic 

213 polynomial fit to the original signal B over windowed time interval t' E [t — A/2, t + A/2], 

214 and T is the time scale not exceeding half of the window length A. The partial derivative 

215 in the above equation is evaluated from the local least-square hnear regression slope of 

216 the Sg{T) dependence in the log - log coordinates for each sliding window. 

217 To our knowledge, the continuous scalogram technique defined by eq. (4) has not been 

218 used in space or turbulence studies before and is introduced in this paper for the first 

219 time. 

220 In this work, we focus on the analysis of magnetic field modulus fluctuations {B = 

221 1 5 (t) I ) providing information on the spectrum of parallel fluctuations of the magnetic field. 

222 These fluctuations are known to be sensitive to ion kinetic effects above the ion spectral 

223 break kpi ~ 1, and represent a distinctive signature of nonlinear compressible cascade 

224 [Alexandrova et al, 2008]. Anisotropic analysis of Mercury's magnetic turbulence, which 

225 will deliver a physically more accurate picture of ion-scale cascades in various Hermean 

226 regions, is left for future research. 

4. Results and Discussion 

4.1. Overview of scaling regimes 

227 Fig. 3 presents the results of the windowed SF analysis of magnetic field fluctuations 

228 observed during the MESSENGER'S first flyby. The studied signal (magnetic field mag- 

229 nitude) is shown on the top panel. The dashed (solid) vertical lines mark the times of 
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230 inbound and outbound crossings of the bow shock (magnetopause) boundaries positioned 

231 according to Slavin et al. [2008]. The dotted vertical hnes show approximate locations 

232 of the outer foreshock boundary identified from our analysis. Upstream of this bound- 

233 ary, the magnetic fluctuations have a quasi-stationary structure of an ambient solar wind 

234 turbulence which is significantly perturbed inside the foreshock region. 

235 Fig. 3(b) shows a stack plot of four time-dependent SF exponents {q — 1 — A, window 

236 size A = 100s) as observed in several selected r channels. During most of the time, the 

237 exponents obey the normal hierarchy with (4^ > d characteristic of a stochastic noise. 

238 There are several noticeable excursions from this rule (marked by arrows) signaling the 

239 presence of transient ULF wave packets as discussed later in the text, the strongest one 

240 being the periodic oscillation in the r channel 0.5-1.0 s detected soon after the closest 

241 approach [Boardsen et al, 2009a, b]. 

242 The bottom panel (Fig. 3(c)) shows the second-order scalogram C2(''', ^) computed for 

243 the same magnetic signal. The local values of the ion crossover scale Tj estimated from the 

244 condition C2 ~ 1 demarking the fluid and ion-kinetic ranges of magnetic turbulence (see 

245 Fig. 1) are plotted with the dashed-dotted line, along with the local proton gyro period 

246 (solid line). The scalogram conflrms the existence of transient ULF wave activity (seen 

247 as pairs of vertically arranged red and blue spots) in several regions visited during the 

248 flyby. More importantly, it shows that the ion crossover scale undergoes a dramatic reor- 

249 ganization during the magnetospheric portion of the flyby, suggesting that even relatively 

250 large-scale plasma motions in this Hermean region should be affected by ion kinetics. To 

251 visuahze this effect, the color coding in Fig. 3(c) is adjusted so that the "kinetic" range 

252 of values of the second-order SF exponent (^2 > 1) is painted in red and the "fluid" range 
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253 (0 < C2 < 1) is in green. The red color clearly prevails inside the magnetospheric cavity. It 

254 can be seen that the interval of scales involved in the kinetic regime grows systematically 

255 as MESSENGER passes through the dusk magnetosheath, and it rapidly expands (by at 

256 least an order of magnitude) during the inbound magnetopause crossing. The outbound 

257 magnetopause crossing is accompanied by an abrupt decrease of r^. The ion crossover 

258 scale remains well above the local proton cyclotron period inside the magnetospheric cav- 

259 ity confirming the presence of strong FLR effects in the Mercury's magnetosphere, in 

260 agreement with numerous previous theoretical predictions (see e.g. Glassmeier and Esp- 

261 ley [2006]; Delcourt et al. [2007]; Blomberg et al. [2007]; Travnicek et al. [2009]; Sundberg 

262 et al. [2010] and refs. therein). 

4.2. Comparative portraits of Hermean plasma structures 

263 Fig. 4 shows the detailed shape results of second-order structure functions describing 

264 the magnetic turbulence in several key plasma regions visited by MESSENGER during 

265 its first flyby. The left (right) columns of plots represent the inbound (outbound) SF 

266 measurements. The discussion below follows the order in which plasma formations were 

267 traversed by the spacecraft. 

268 SoIelt wind on both disk and dawn sides of Mercury exhibits classical signatures of 

269 large-scale fluid cascade coexisting with kinetic-scale turbulence. Some variability of low- 

270 frequency SF exponents seen in Fig. 3(b) can be due to the inherent intermittency of the 

271 solar wind flow [Roberts et al, 1992; Borovsky, 2010]. The structure functions have a 

272 crossover at r fa 0.5-1.0 s (Fig. 4(a-b), black curves). The exponent above this scale is 

273 reasonably close to the Kolmogorov's law (/3 = 5/3, ^2 = 2/3) in the outbound solar wind, 

274 and is more consistent with the Iroshnikov - Kraichnan scaling ansatz (/? = 3/2, (^2 — 1/2) 
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275 [Biskamp, 2003] in the outbound solar wind. For r < 0.5s, the SF slopes are considerably 

276 steeper. The value C2 > 1 observed in this range of scales is indicative of the ion-kinetic 

277 regime, and it implies that the power spectral density of the magnetic fluctuations scales 

278 as k"^ with /? = C2 + 1 > 2. 

279 Compared to the inbound solar wind region (Fig. 4 (a)), the outbound solar wind mea- 

280 surements (panel (b)) are less stable, and they exhibit a more pronounced ion- kinetic 

281 component propagating toward larger r , possibly reflecting wave turbulence initiated in 

282 the foreshock region upstream of the bow shock. 

283 The foreshock region contains a strongly inhomogeneous turbulent environment filled 

284 with transient packets of quasi-periodic oscillations and high-frequency stochastic noise. 

285 During the inbound portion of the flyby, the solar wind structure undergoes an abrupt 

286 change at the upstream foreshock boundary which first affects the kinetic scales of the 

287 turbulent spectrum leaving the larger (MHD) scales almost unperturbed, see dashed red 

288 curve in Fig. 4 (a). After this magnetic fluctuations reorganize themselves across the entire 

289 T range (sohd red curve, same Figure). The repetitive decreases of short-scale SF expo- 

290 nents (marked with arrows in Fig. 3(b)) indicate that the spacecraft has flown through 

291 several regions of ULF wave activity in both dusk- and dawn-side foreshocks as discussed 

292 above. 

293 Our observations show that the Mercury's coUisionless foreshock possesses a well- 

294 developed macrostructure possibly associated with ULF waves and discontinuities gen- 

295 erated by backstreaming ions [Fairfield, 1991; Omidi et ai, 2006]. Some of the detected 

296 intermittent structures can be due to hot flow anomalies upstream of the bow shock such 

297 as the ones found recently at Venus [Slavin et ai, 2009b]. Three-dimensional kinetic 
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298 simulations predict that the outbound foreshock may contain beams of plasma directed 

299 from Mercury's bow shock back upstream against the solar wind flow, resulting in a com- 

300 plex regime of wave-particle energy exchange manifested in long- wavelength beam-driven 

301 oscillations [Travnicek et al, 2009]. 

302 The magnetosheath is dominated by intermittent kinetic fluctuations with nonlinear 

303 spectrum (not shown) converging to the ion-kinetic regime for r < 1 s, reminis- 

304 cent of the turbulence in the terrestrial magnetosheath as observed by Cluster spacecraft 

305 [Yordanova et al, 2008]. The stochastic component is mixed with transient episodes of 

306 ULF oscillations of various frequencies, the most intense ULF episode being observed in 

307 the inbound magnetosheath during 18:20-18:22 UT (shown by arrow in Fig. 3(b)) at 

308 a time scale of about one third of the local proton gyroperiod. The reorganization of 

309 magnetic fluctuations at the magnetosheath entry has begun from large scales (dashed 

310 blue curve in Fig. 4(c)) and involved shorter scales in about 15 minute after the inbound 

311 bow shock crossing. A fully-developed broad-band kinetic turbulence obtains in the near- 

312 magnetopause region of enhanced rms variability (solid blue curve, same panel). A similar 

313 sequence of events (in the reversed order) was observed while crossing the outbound bow 

314 shock (Fig. 4(d)) which also contains intense packets of ULF oscillations (Fig. 3 (b-c)). 

315 Hybrid simulations show that downstream of the bow shock. Mercury's plasma is 

316 marginally stable with respect to mirror and cyclotron instabilities producing large- 

317 amplitude compressible waves [Travnicek et al., 2009]. The same study suggests that 

318 the outbound magnetosheath can be also prone to fire-hose instabilities. It remains to be 

319 verified whether the ULF episodes present in our results for the Hermean magnetosheath 

320 are associated with some of these instability mechanisms. 
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321 During inbound magnetosheath observations, one flux transfer event (FTE) at 

322 UT=18:36:21-18:36:25 lias been documented [Slavin et al, 2008, 2010a]. PTEs in tfie 

323 magnetoslieatli are produced by localized magnetic reconnection between the interplane- 

324 tary and planetary magnetic fields, and are seen as passages of helical magnetic structures 

325 with a characteristic bipolar By signature encompassing a core region of an increased field 

326 magnitude. The magnetic fluctuations during the FTE (yellow curve in Fig. 2(c)) differ 

327 from those characterizing average conditions in the surrounding part of Mercury's mag- 

328 netosheath. Based on the shape of the second-order SF showing a nearly Kolmogorovian 

329 scahng, the FTE has launched a partly-developed fluid cascade modulated by a quasi- 

330 periodic distortion at r ~4-6 s consistent with the time scale of the bipolar By signature 

331 [Slavin et al., 2008]. Our analysis also hints at the possibility of multiple FTEs and/or 

332 thin current sheets in the inbound magnetosheath during 18:24 - 18:37 UT (outlined by 

333 rectangle in Fig. 2(b)). Their presence is suggested by the anti-correlation between the 

334 SF exponents measured at r = 5-15 s and r = 2.5-5 s implying several transient features 

335 on a ~5 second time scale. 

336 Mercury's magnetosphere reveals a rich diversity of scahng regimes most of which 

337 are shaped by kinetic-type fluctuations. The inbound magnetopause crossing is marked 

338 by a rapid transition from the magnetosheath turbulence characterized by relatively nar- 

339 row range of kinetic behavior, to a developed kinetic turbulence described by C2 = 3/2 

340 (indicative of ion-kinetic cascade [Schekochihin et al, 2007]) over broad range of scales. 

341 The rather high upper time scale limit of ion-kinetic turbulence in the Kelvin- Helmholz 

342 instability region (Fig.4(2), dashed green line) matches the average period (~ 20s) of vor- 
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343 tex rotations [Slavin et ai, 2008] and therefore does not necessarily represent an intrinsic 

344 fluid crossover such as the one observed in the solar wind (Fig. 4(a)). 

345 The equatorial plasma sheet (same panel, solid green line) displays ion-kinetic turbulent 

346 scaling across the entire studied range of r. This is quite different from the behavior of the 

347 terrestrial current sheet outside the reconnection region. In the geotail, the dissipation 

348 and kinetic effects usually play a leading role at r < 1 s while larger scales tend to be 

349 controlled by an intermittent fluid cascade with /3 < 2.5 (^2 < 1-5), see Voros et al. [2006] 

350 for a brief review. One can infer that the dynamics of the central plasma sheet in the 

351 Hermean magnetosphere is strongly affected by non-MHD effects introduced by finite sizes 

352 of cyclotron orbits of the constituting ion species, in agreement with earlier theoretical 

353 predictions (see e.g. Delcourt et al. [2007]). As discussed in section 4.3, CCS turbulence 

354 is consistent broadly with a very quiet, thick plasma sheet reported by Slavin et al. [2008] 

355 based upon the large magnetic field in the depressed equatorial tail and the northward 

356 IMF Bz in the solar wind during the first flyby. The key open problem is the nature of 

357 the waves in the dissipation range - whether the turbulent energy is deposited in the form 

358 of kinetic Alfven waves or whistler waves [Eastwood et al., 2009]. Without simultaneous 

359 electric and magnetic field measurements, this question may not have a definite answer. 

360 The fiuid component of the magnetospheric turbulence can be reliably identified only 

361 during the near-Mercury portion of the flyby, namely during the flrst diamagnetic decrease 

362 encountered in the inner magnetosphere [Slavin et al., 2008], gray line in Fig. 4(f). Based 

363 on the analysis of a similar region in the terrestrial magnetosphere [Uritsky et al., 2010b; 

364 Liu et al., 2011; Panov et al., 2010], these fluctuations can manifest transient velocity 

365 and magnetic fleld shears due to reconnection - driven sunward flow bursts in the plasma 
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366 sheet. The flows are expected to stir turbulent vortices at the inner edge of the plasma 

367 sheet where the sunward convecting plasma sheet ions encounter the stronger planetary 

368 dipole magnetic field and are quickly decelerated [Shiokawa et ai, 1998]. At smaller radial 

369 distances in the Earth's magnetotail, fluid turbulence is suppressed due to a stabilizing 

370 effect of the dipole magnetic field [Stepanova et al., 2009, 2011]. A much weaker dipole 

371 field at Mercury apparently allows turbulent vortices to penetrate closer to the planetary 

372 surface. 

373 Physical interpretation of the ULF wave activity observed after the closest approach 

374 (see Fig. 3) remains a challenging task. Although the frequency of these waves is close to 

375 the local proton gyro frequency, their mixed polarization, with a large amount of right- 

376 hand polarized packets [Boardsen et al., 2009a, b], does not fit the conventional picture 

377 of ion-cyclotron resonant instability. A series of higher harmonics detected in the inner 

378 magnetosphere during the first flyby suggests that the observed ULF oscillations may in 

379 fact represent magnetosonic waves driven locally by a non-maxwellian proton distribution 

380 [Anderson et al, 2011]. 

381 The last (third) diamagnetic decrease traversed by MESSENGER before the outbound 

382 magnetopause crossing demonstrates no signatures of fluid or MHD scaling (Fig. 4(f), 

383 dashed orange line) . The shape of the structure function of magnetic field modulus fluctu- 

384 ations in this region is close to that in the first diamagnetic decrease. A more sophisticated 

385 anisotropic analysis (to be published elsewhere) shows a distinct similarity between ra- 

386 dial component of magnetic fluctuations in this region and at the adjacent magnetopause 

387 boundary (solid blue line in Fig. 4(i)). This observation provides an indirect support to 

388 the hypothesis by Slavin et al. [2008] who classified the third diamagnetic decrease as an 
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389 ion boundary layer compatible with the gyro-radius of sodium pickup ions accelerated in 

390 the magnetosheath, and described this layer, together with the outbound magnetopause, 

391 as an integrated double-magnetopause structure. On the other hand, simulations show 

392 that temperature anisotropy in this region can be regulated by proton mirror and proton 

393 cyclotron instabilities [Travnicek et ai, 2009]. Judging from the spectral amplitude of 

394 magnetic fluctuations in the vicinity of proton gyro frequency [Anderson et al, 2011], 

395 kinetic effects in the IBL are likely to have much higher growth rates and/or saturation 

396 levels than those in the inner magnetosphere. A proper kinetic treatment of the IBL region 

397 involving a multi-ion plasma composition and accurate resolution of relevant instability 

398 scales seems to be necessary for understanding the underlying physics of this complex 

399 plasma structure. 

400 The large-scale behavior of the SF plots presented in Fig. 4 enables an indirect verifica- 

401 tion of the stationarity of the studied data segments. Following the approach proposed by 

402 Matthaeus and Goldstein [1982], the stationarity of magnetic fluctuations can be tested 

403 based on the ergodic theorem for stationary random processes [Monin and Yaglom, 1975]. 

404 In its simplest version, the theorem states that the time average of B obtained over 

405 subintervals of a limited duration converges to the ensemble average as the length of the 

406 subintervals significantly exceeds the correlation time of the signal. This condition en- 

407 sures so-called weak stationarity of magnetic turbulence, and it tends to be fulfilled in 

408 the interplanetary medium but not necessarily for planetary magnetic fluctuations. The 

409 ergodic convergence poses a restriction on the asymptotic shape of the two-time autocor- 

410 relation function R{t) which must decay as or faster (see Matthaeus and Goldstein 

411 [1982] for details). It can be easily shown that this requirement is violated for ^2 > 1 
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412 since in this case R{t) oc r~" with a = 2 — ^2 < 1 [Carreras et al, 1999; Li, 2010], but 

413 is met for C2 < 1- The large-scale log-log slopes of all the plots in Fig. 4 are smaller 

414 than 1, and so the studied signals are stationarity at least in the weak sense. For some 

415 other locations in Mercury's magnetosphere, however, the situation is not as clear. For 

416 instance, during the last two minutes (18:48-19:00) prior to entering the near- Mercury DD 

417 region, the crossover scale defined by the condition C2 ~ 1 was quite close to the upper 

418 measured time scale, see Fig. 3(c). This and similar regions require a more accurate 

419 stationarity analysis addressing convergence of higher statistical moments, and are not 

420 used for quantitative calculations in the Section 4.3 below. 

421 For comparison purposes. Fig. 5 presents Fourier power spectra of magnetic turbulence 

422 in some of the magnetospheric regions discussed above. As expected, the shape of the 

423 spectra is consistent with the ion kinetic regime 11 (see Fig. 1) described by /3 ~ 2.5, but 

424 is statistically less stable than the shape of the SFs in the same range of scales (Fig. 4). 

425 The power spectra do not resolve the fine low-frequency structure of the studied signals 

426 which is clearly seen in the more robust structure function statistics. Anderson et al. 

427 [2011] have obtained Fourier power spectra of linearly detrended magnetic field data in 

428 the inner magnetosphere and in the IBL region. The spectral densities reported in their 

429 work undergo a rather steep decay characterized by ^ fa 2.5 — 3.0 above the local proton 

430 cyclotron frequency, which is roughly consistent with the shape of the spectra in Fig. 5. 

4.3. Quantitative estimates 

431 Table 2 summarizes the results of our calculations of ion gyro radii and temperatures 

432 using eq.(2-3). The proton and sodium temperatures have been estimated assuming that 

433 the crossover sale is controlled by the gyromotion of protons and Na"*" ions, correspond- 
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434 ingly. We have considered only those plasma regions for which the Taylor frozen-in flow 

435 condition is roughly satisfied and so the linear mapping between the spatial and temporal 

436 domains of analysis is possible. As input parameters, we used the predicted values of bulk 

437 fiow velocities for Hermean plasma environment [Slavin et al., 2008; Baker et al., 2009], 

438 the average flow velocity in a quiet Earth's central plasma sheet [Angelopoulos et al., 1993; 

439 Baumjohann et al, 1989], as well as the values of the proton gyro period Tcp and the ion 

440 crossover Tj obtained from our analysis. 

441 Solar wind parameters on the dusk (SWl) and dawn (SW2) sides of Mercury's magne- 

442 tosphere are roughly consistent with the results of global solar wind simulations [Baker 

443 et al, 2009, 2011] which predict pi ^ 18 km and Tp ^ 10 eV for the flrst MESSEN- 

444 GER's flyby. If we accept that Tj in the solar wind is controlled by the ion inertial length 

445 Aj = c/ujpi, where c is the speed of light and ujpi is the ion plasma frequency, rather than 

446 by Pi (see Sahraoui et al. [2009] and Sahraoui et al. [2010] for details), the turbulence- 

447 based estimates become closer to the simulated values. For a plasma beta (the ratio of 

448 the particle pressure to magnetic pressure) of the order 1, this assumption yields pi ~ 20 

449 (30) km, Tp w 11 (25) eV, and the proton number density rip ~ 65 (40) cm~^ at the dusk 

450 (dawn) flanks. The inbound numbers are in a good agreement with Baker et al. [2009] 

451 as well as with typical solar wind parameters at Mercury orbit [Blomberg et al., 2007]. 

452 The outbound estimates reveal somewhat hotter plasma environment, possibly due to an 

453 extended quasi-parallel foreshock system existing at the dawn side. 

454 The magnetosheath plasma estimates vary greatly with time and position as can be 

455 expected from the behavior of nonstationary Q exponents and the SF scalogram con- 

456 structed for this region (Fig. 3). Sundberg et al. [2010] have reported a characteristic 
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457 proton temperature in the Mercury's magnetosheath of about 700 eV obtained by scaling 

458 the terrestrial magnetosheath proton temperature with the scaling factor given in Slavin 

459 and Holzer [1981]. Our temperature estimates for the inbound (MS2) and outbound 

460 (MS3) magnetosheath regions adjacent to the magnetopause are significantly below this 

461 predicted value. However, the proton temperature in many other magnetosheath locations 

462 exceeds the prediction. For example, the flux transfer event in the dusk magnetosheath 

463 is characterized by Tp 1200 eV and the ion Larmor radius of ~ 240 km, or about 20% 

464 of the size of the FTE as estimated by Slavin et al. [2008]. The flux rope topology of this 

465 FTE should therefore be considerably affected by FLR effects. 

466 Hybrid simulations of the first fiyby [Travnicek et al, 2009] demonstrate a steep in- 

467 crease (by a factor of 3) in the ion temperature during the inbound magnetopause cross- 

468 ing, accompanied by a noticeable drop in the plasma density. They also predict strong 

469 temperature gradient at the dawn magnetopause. Both transitions are clearly present in 

470 our results showing a sharp increase (decrease) of the estimated Tj and pi values during 

471 the inbound (outbound) magnetopause crossings. These transitions are also captured by 

472 the dramatic growth of at the dusk boundary and the decay of the this parameter on 

473 the dawn side (see Fig.3(c)). 

474 The conditions at the magnetopause boundary suggest a significant contribution from 

475 FLR effects, with the largest pi ~ 0.1i?M observed just inside the inbound magnetopause. 

476 Due to final gyro orbits, the dusk side magnetopause can be either less stable than the 

477 dawn magnetopause, or not stable at all [Glassmeier and Espley, 2006]. The pronounced 

478 signatures of KH activity observed at the inbound magnetopause during the first fiyby 

479 confirmed the prediction [Slavin et al., 2008; Sundberg et al., 2010]. The possibility of 
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480 Kelvin - Helmholtz vortices on the opposite side of Mercury's magnetopause remains 

481 controversial. Using an FLR extension of the ideal MHD, Glassmeier and Espley [2006] 

482 have shown that the smallest KH- unstable wavenumber A^j„ = 8TTr]/\6v\ at the dawn 

483 magnetopause is controlled by the the kinematic viscosity rj — p^uJci/4: and the magnitude 

484 \Sv\ of the velocity shear. Our measurements suggest that both the dusk and the dawn 

485 flanks of Mercury's magnetosphere can be prone to KH instability. By plugging the 

486 plasma parameters of dawn magnetopause (Table 2) into the above expressions, we obtain 

487 rj ~ 10^ — 10^ m^/s and Xmin of the order of 100 km. This implies that the KH growth 

488 rate in this region can be positive for a wide range of wavenumbers, in agreement with 

489 recent theoretical results [Sundherg et al, 2010]. 

490 The narrow-band ULF wave packets observed by MESSENGER between the closest 

491 approach and the outbound magnetopause [Slavin et al, 2008] should be also strongly 

492 affected by finite gyro radii. The frequency of these waves has been found to be close 

493 to the local He"*" cyclotron frequency [Boardsen et al., 2009a, b], which is by an order of 

494 magnitude larger than the frequency corresponding to the ion crossover scale in this region 

495 (tj ~ 10 s). The generation mechanism of these Hermean wave packets is very likely to be 

496 kinetic. The measured Tj value is comparable with local cyclotron period of sodium ions 

497 (~ 20 s) hinting at their involvement in the cross-scale coupling processes in the studied 

498 magnetospheric region. This is quite different form Earth's magnetosphere where ULF 

499 waves tend to have frequencies well below all relevant gyro frequencies [Blomberg et al, 

500 2007]. 

501 According to our investigation, the cross-tail current sheet (CCS) plasma population at 

502 Mercury can be significantly denser and cooler than the one typically observed at Earth. 
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503 The size of the gyro radius reported in Table 2 is also by a factor of two smaller than the 

504 one inferred from the measurements performed by MESSENGER Fast Imaging Plasma 

505 Spectrometer (FIPS) {Tp ~ 2 x 10^ K, or 170 eV, according to Raines et al. [2010]) which 

506 yield pi ~ 120 km. Pressure-balance arguments [Slavin et al., 2011] provide a current 

507 sheet plasma beta of ~ 5 for the studied flyby. Using this beta value, we obtain Up ~ 40 

508 cm~^, which exceeds FIPS estimates [up ~ 1-10 cm~^). 

509 Our CCS results are not completely unexpected considering the northward IMF orien- 

510 tation during the studied time interval. A super-dense, cool plasma sheet material similar 

511 to the one reported here has been sighted in the terrestrial magnetosphere during extreme 

512 geomagnetic calm intervals characterized by steady northward IMF (see Borovsky and 

513 Steinberg [2006] and refs therein). Such calm intervals may be important for precondi- 

514 tioning the magnetosphere for subsequent geomagnetic perturbations. Cool plasma can 

515 be an effective contributor to the inner magnetosphere since cool plasma sheet particles 

516 are less subject to gradient curvature drift and can be convected deeper into the dipole re- 

517 gion and therefore producing greater adiabatic pressure increase compared to hot particles 

518 [Borovsky and Steinberg, 2006]. 

519 By applying the simple 1/r^ correction (where r is the distance from the Sun, see 

520 Ogilvie et al. [1977]) to the proton densities of ~ 3.0 — 4.5 cm~^ reported by Borovsky 

521 and Steinberg [2006] for the dense terrestrial plasma sheet, we expect the corresponding 

522 Hermean values to be in the range rip ~ 30 — 50 cm~^, which agrees with the CCS density 

523 measured in our study {up ~ 40 cm~^). Mukai et al. [2004] have extrapolated data-driven 

524 terrestrial plasma sheet models by Terasawa et al. [1997] and Tsyganenko and Mukai 

525 [2003] to account for the substantially different interplanetary environment at Mercury 
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526 and get a baseline for the plasma analyzer onboard the upcoming BepiColombo mission. 

527 Using this empirical extrapolation, the solar wind density ~ 65 cm~^ translates into the 

528 plasma sheet density of ~ 1 — 20 cm~^ (see Fig. 2(b) of Mukai et al. [2004]). The upper 

529 limit of 20 cm~^ corresponds to the cold and dense state of the Hermean plasma sheet 

530 and is comparable with our density estimate. 

531 If the estimates provided in Table 2 are correct, the relatively small ion scales in the 

532 Mercury's CCS could help explain the short characteristic substorm time scale of ~ 1 — 3 

533 min [Baumjohann et al, 2006; Slavin et al, 2009c, 2011] on this planet. For Hermean 

534 substorms to be this short, tail reconnection at Mercury has to be extremely fast and 

535 intense [Blomberg et a/., 2007]. The reconnection rate is largely controlled by the current 

536 sheet thickness which is of the order of the ion skin depth \ ~ 35 km, based on our 

537 assessment. For the convective inflow speeds of several hundreds kilometers per second, 

538 the transition time of this depth would be a few 100 milliseconds. Furthermore, if the 

539 reconnection at Mercury proceeds inside the ion diffusion region on electron inertial scale 

540 which we estimate to be ~ 1 km, the transition time could be as little as 10 ms, making 

541 fast impulsive reconnection possible and perhaps inevitable. 

542 In the absence of sufficiently intense tail lobe loading, no actual substorm activity was 

543 observed during the studied fly by. In agreement with this fact, our analysis suggests 

544 that proton trajectories in Mercury's current sheet were nearly adiabatic. The effects of 

545 magnetic moment scattering in thin current sheets can be conveniently measured by the 

546 adiabaticity parameter k, introduced by Buchner and Zelenyi [1989]. By definition, k is 

547 the square root of the ratio of the smallest field-of-line curvature radius to the largest ion 

548 Larmor radius. For particles traveling through a field reversal, the condition «; > 3 ensures 
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549 adiabatic behavior; for k < 3, the magnetic moment scattering is responsible for particle 

550 injection into the loss cone [Sergeev et al, 1983], with a possibility of parametric "islands" 

551 of quasi-adiabatic behavior at very small k values [Delcourt et ai, 2006]. Assuming that 

552 the field line curvature radius is of the order of 1 Rm — 2440 km and that the measured 

553 Pi approximates the maximum relevant Larmor radius, we get k ~ 7. The result is well 

554 above the transitional value k — 3, and is a signature of adiabaticity. To unfreeze the 

555 magnetic flux and initiate reconnection, a much more stretched magnetotail configuration 

556 would be required. Such a configuration has been reached during the second and the third 

557 messenger's flybys which revealed a rather strong dayside and nightside reconnection 

558 activity accompanied by intense loading and unloading events [Slavin et al., 2009c, 2010b]. 

559 One more indication of a stable state of the current sheet is its relatively low Reynolds 

560 number Re evaluated using {Ref/^ r^Lji [Warhaft, 2002], where L and £ are the largest 

561 and the smallest scales of the inertial range cascade, respectively. In the CCS case, L is 

562 defined by the size of the flow channel of the radial convective plasma transport. In the 

563 terrestrial plasma sheet, L is about 10% of the width of the plasma sheet, or ~ 2 Earth 

564 radii [Nakamura et al., 2004], £ ft! 50 km, and so Re ~ 1600 [Voros et al., 2006], which is 

565 indicative of a marginally stable regime at the edge of turbulence. At Mercury, due to a 

566 smaller planetary size, the estimated Reynolds number is much lower. If we assume that 

567 the BBF channel width at Mercury is also 10% of the width of the tail, then L = 0.5Rm 

568 in the Hermean magnetosphere. Alternatively, using the scaling factor of 8 given by the 

569 ratio of terrestrial to Hermean magnetospheric sizes in units of the respective planetary 

570 radii [Ogilvie et al, 1977], one can argue that the flow channel of 2 Earth radii becomes 

571 ~ 0.25i?M at Mercury. By using L = 0.25 — 0.5Rm and substituting pi = 60 km from 
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572 Table 2 as a proxi to we arrive at Re ~ 22 — 55. This fairly small Reynolds number 

573 suggests a predominantly laminar regime of plasma flow in the Hermean cross-tail current 

574 sheet, consistent with the shape of the SF in this region (Fig. 4(e)) which reveals a rather 

575 limited interval of scales of fluid cascade, if any at all. 

5. Conclusion 

576 We have presented the results of a first investigation of magnetic fluctuations in the 

577 near-Mercury space environment. Our main findings can be summarized as follows: 

578 (1) Turbulent conditions in the solar wind during the studied fly by were close to stan- 

579 dard, with well-developed MHD and ion-kinetic components, consistent with the results 

580 reported by Korth et al. [2010] for MESSENGER'S solar wind observations. 

581 (2) Foreshock plasma at Mercury is populated with transient oscillatory perturbations 

582 organized over large spatial distances. This macrostructure can be associated with magne- 

583 tosonic waves and ULF plasma modes generated by field-aligned ion beams as proposed 

584 earlier [Omidi et al, 2006]. The low frequency fluctuations in the foreshock have no 

585 obvious association with any known type turbulent cascade. 

586 (3) The magnetosheath turbulence is dominated by intermittent kinetic-scale fluctu- 

587 ations, in agreement with similar observations at Earth. Judging from a single FTE 

588 observation in the magnetosheath, traveling flux ropes can be a source of enhanced low- 

589 frequency turbulence in this plasma region. 

590 (4) Turbulence in Mercury's magnetosphere is strongly influenced by finite gyroradius 

591 effects, with fluid-type energy cascades playing secondary or no part in most of the re- 

592 gions inside the magnetospheric cavity, which supports earlier theoretical predictions and 

593 simulation results. 
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594 (5) Stochastic properties of the central current sheet in Hermean magnetotail speak 

595 in favor of its relatively stable global configuration consistent with the steady northward 

596 IMF driving and the absence of noticeable substorm activity during the first flyby. 

597 Overall, our results show, for the first time, that turbulence in the Hermean magneto- 

598 sphere as well as in the surrounding space region is strongly affected by non-MHD effects 

599 introduced by finite sizes of cyclotron orbits of the constituting ion species. We conclude 

600 that kinetic effects may play a critically important role in the Mercury's magnetosphere up 

601 to the largest resolvable time scale (~ 20 s) imposed by signal nonstationarity. However, 

602 the prevalence of turbulence signatures of kinetic processes does not necessarily mean that 

603 the latter are determining the structure of Mercury' magnetic field. Rather, our results 

604 indicate that these kinetic processes need to be identified, and their potential influence on 

605 Hermean magnetosphere need to be understood. A more sophisticated statistical analysis 

606 addressing multiscale anisotropic properties of magnetic turbulence formed under differ- 

607 ent solar wind driving conditions will be required to clarify physical mechanisms of these 

608 effects and their influence on other Hermean processes such as e.g. a tail reconnection, 

609 plasma transport, generation of field-aligned currents, and ULF wave activity. These and 

610 related tasks outline a fruitful field of future research. 
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Power spectrum Structure function 




Figure 1. Schematic diagram sliowing the typical shape of the power spectrum and the 
structure function of magnetic turbulence involving abrupt changes (scaling crossovers) 
of the power-law exponents /3 and ( at the ion and electron gyro scales. The hierarchy 
of the exponents in the fluid, ion- and electron-kinetic ranges of scales (regions I, II and 
III, correspondingly) is shown. The ion crossover characterized by kpi ~ 1 is the main 
subject of this work. 
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MESSENGER Mercury Encounter 1 14 January, 2008 
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Figure 2. MESSENGER'S first flyby trajectory overlapped with the average bow shock 
and magnetopause boundaries obtained using crossing information from five available 
fiybys [Slavin et al, 2009a]. Color bars show some of the studied regions, see Table 1 for 
full description. 
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Table 1. Intervals of analysis used for computing structure functions in Fig. 4-5. 



Notation Time Description 



SWl 


17:10:00- 


■17:40:00 


Unperturbed solar wind at the dusk side 


FSl 


17:45:00- 


•17:48:00 


Outermost dusk-side foreshock region 


FS2 


18:05:30- 


■18:08:30 


Innermost foreshock region near the dusk-side bow shock 


MSI 


18:09:00- 


■18:12:00 


Outermost magnetosheath at the dusk flank 


PTE 


18:36:00- 


■18:37:00 


One-minute interval involving flux transfer event 


MS2 


18:39:00- 


■18:42:00 


Innermost magnetosheath contacting the inbound magnetopause 


KH 


18:43:00- 


■18:46:00 


Kelvin - Helmholtz vortices at the dusk magnetopause 


CCS 


18:47:00- 


■18:49:00 


Cross-tail current sheet 


DD 


19:00:00- 


■19:03:00 


First diamagnetic decrease encountered in the inner magnetosphere 


IBL 


19:11:00- 


■19:14:00 


Ion boundary layer adjacent to the outbound magnetopause 


MS3 


19:14:30- 


■19:17:00 


Innermost magnetosheath observed after exiting the magnetosphere 


MS4 


19:17:00- 


■19:18:30 


Outermost magnetosheath adjacent to the outbound bow shock 


FS3 


19:19:30- 


■19:22:30 


Innermost outbound foreshock region 


FS4 


19:42:00- 


■19:45:00 


Outbound foreshock adjacent to the unperturbed solar wind 


SW2 


19:52:00- 


■20:22:00 


Unperturbed solar wind at the dawn side 
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Figure 3. Windowed higher-order SF analysis of magnetic field fluctuations recorded during the first 
MESSENGER flyby. (a) Time series of the total field magnitude B. (b) Time-dependent structure function 
exponents (C^\ - black, <^2 - blue, (^3 - green, C,a - red) estimated within five ranges of temporal scales of B 
variability. Black arrows show episodes of ULF activity discussed in the text, (c) Continuous second-order SF 
scalograms C,2(j,t) computed for the same signal. The red color corresponds to the fully developed ion-kinetic 
turbulent cascade with /3 = 2.5 and Qi = 1.5. Black solid line overplotted with the scalogram shows the local 
proton cyclotron period; dashed-dotted curve is an approximate ion crossover time scale evaluated from on the 
^2 ~ 2 condition using the SF analysis. The vertical solid, dashed, and dotted lines mark the inbound and 
outbound positions of the magnetopause, bow shock, and foreshock, correspondingly. 
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Figure 4. Second-order structure functions of magnetic field modulus fluctuations 
characterizing MESSENGER'S crossings of key Hermean plasma boundaries and struc- 
tures. Left and right panels show inbound and outbound encounters, correspondingly. 
Proton cyclotron periods are shown with vertical lines of matching color and pattern on 
the bottom of each panel. Tilted straight lines representing theoretical slopes for the fully 
developed fluid (K41, C2 = 2/3) and Hall MHD (HMHD, C2 = 3/2) scaling regimes are 
added for reference to each plot. Notations and time limits of the studied regions are 
explained in Table. 1. 
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Figure 5. Fourier power spectra of several magnetospheric regions showing broad-band 
kinetic fluctuations consistent with Hall MHD interpretation. The spectral power law 
index (3 is related with the second-order SF exponent as /3 = C2 + 1- 
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Table 2. Estimated plasma parameters in selected regions of Hcrmcan magnctosphcre: 

Tcp - proton cyclotron period; Vq - typical bulk fluid velocity; Tj - ion crossover scale 

corresponding to the transition between the fluid- and kinetic-like behavior of the structure 

function; pi - ion gyroradius obtained using eq.(2); Tp (Tno) ~ temperature of protons 

(Na"*" ions) evaluated from (3); Rm ~ 2440 km - radius of Mercury. See Table 1 for region 

notations. 

Region Tcp, s vo, km/s n, s pi, km pi/Rn Tp, eV T^a, eV 



SWl 3.9 

PTE 4.6 

MS2 5.5 

KH 5.9 

CCS 6.0 

IBL 1.0 

MS3 1.6 

SW2 3.2 



450 0.4 

300 5.0 

300 1.5 

150 10.0 

50 7.0 

150 3.0 

300 0.4 

450 0.5 



30 0.01 

240 0.10 

70 0.03 

240 0.10 

60 0.03 

70 0.03 

20 0.01 

35 0.01 



25 

1200 50 

75 3 

700 30 

40 2 

2300 100 

60 3 

55 
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